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Abstract 


Tests were conducted in a Mach 6 flow to examine the reattachment process of an axisymmetric 
free shear layer associated with the near wake of a 70 deg half angle , spherically blunted cone with a 
cylindrical after body . Model angle of incidence was fixed at 0 deg and free-stream Reynolds numbers 
based on body diameter ranged from 0.5xl0 6 to 4xl0 6 . The sensitivity of wake shear layer transition 
on reattachment heating was investigated . The present perfect gas study was designed to compliment 
results obtained previously in facilities capable of producing real gas effects . The instrumented blunted 
cone model was designed primarily for testing in high enthalpy hypervelocity shock tunnels in both this 
country and abroad but was amenable for testing in conventional hypersonic blowdown wind tunnels 
as well Surface heating rates were inferred from temperature - time histories from coaxial surface 
thermocouples on the model forebody and thin film resistance gages along the model base and 
cylindrical after body . General flow feature (bow shock , wake shear layer, and recompression shock) 
locations were visually identified by schlieren photography. Mean shear layer position and growth were 
determined from intrusive pitot pressure surveys. In addition , wake surveys with a constant 
temperature hot-wire anemometer were utilized to qualitatively characterize the state of the shear layer 
prior to reattachment. Experimental results were compared to laminar perfect gas predictions provided 
by a 3-D Navier Stokes code (NSHYP). 

Shear layer impingement on the instrumented cylindrical after body resulted in a localized 
heating maximum that was 21 to 29 percent of the forebody stagnation point heating. Peak heating 
resulting from the reattaching shear layer was found to be a factor of 2 higher than laminar 
predictions, which suggested a transitional shear layer. Schlieren flow visualization and fluctuating 
voltage time histories and spectra from the hot wire surveys across the shear layer substantiate this 
observation. The sensitivity of surface heating to forebody roughness was characterized for a 
reattaching shear layer. For example, at Reoo t d = 4x10 6 , when the shear layer was transitional, the 
magnitude of peak heating from shear layer impingement was reduced by approximately 24 percent 
when transition grit was applied to the forebody. The spatial location of the local peak, however, 
remained unchanged. 


Nomenclature 

Ch Stanton number, 

CfJ [rooUoo(rht -h^)] 
d diameter, (in.) 

d boundary layer thickness, mm 

h enthalpy, (J/kg) 

k nominal roughness diameter, (mm) 

M x free stream Mach number 

pt,i tunnel stagnation pressure, N/m 2 

Pt,2 free stream pitot pressure, N/m 2 


r wall recovery factor 

R radius, (in.) 

S distance along model surface 

t time, (sec) 

Reoo unit free stream Reynolds number, 

(m 1 ) 

T temperature, (K) 

u velocity, (m/sec) 

x axial distance from model base 
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Subscripts 

00 free-stream conditions 

j juncture 

m model forebody 

n model nose 

ref laminar reference heating 

s sting 

stag reference heating 

t reservoir stagnation conditions 

w wall 


Introduction 

NASA has examined the concept of 
aerobraking to satisfy planetary mission 
requirements. Aerobraking as defined by 
Jones, 1987, involves the use of aerodynamic 
forces acting on a spacecraft surface to 
decelerate. One proposed application of this 
technique for a planetary mission would be the 
placement of a science payload (Tauber et al, 
1993; Mitcheltree, 1994) on or into the surface 
following a direct entry from an interplanetary 
flight. 

One concept of this type, shown in Fig. 1 , 
envisioned the use of a blunt umbrella shaped 
forebody to provide the aerobraking surface. 
The spacecraft design for a planetary direct 
entry approach can be similar. In this mission 
scenario, the desired payload and/or 
instrumentation package to be delivered to the 
surface would be placed behind the aerobrake 
to protect it from the intense heat generated 
during atmospheric entry. The Mars 
Pathfinder and Mars 98 missions under the 
NASA Discovery Program are examples of 
future planned missions utilizing this type of 
technology. 

The accurate prediction of the near wake 
flow structure associated with blunt aerobrake 
concepts is an essential part of the design 
process as payload constraints are often 
imposed by wake closure. Proper positioning of 
the payload will be critical to insuring 
aerodynamic stability of the spacecraft during 
entry and to avoid thermal damage from 
localized near wake phenomenon. Although it 
is generally recognized that convective heating 
rates in the base region are low (Gnoffo, 1992), 
localized maxima can occur if the free shear 
layer that separates from the corner of the 


aerobraking surface (forebody) impinges on the 
after body (Gnoffo, 1992; Wells, 1990; Dye, 
1993). The successful design of an aerobraking 
vehicle therefore requires an understanding of 
the physics associated with hypersonic, blunt- 
body, near-wake flowfields. 

Understanding the physical aspects of the 
compressible free shear layer is of general 
interest as it can also be found in other 
practical situations. For example, internal 
flows involving axial flow combustors (King et. 
al., 1990; Herrin, 1995) often result in wake 
shear layers that are produced by fuel injectors. 
The growth and stability of these shear layers 
has a direct impact on the fuel-air mixing 
efficiency in scramjet combustors (e.g. 
Macaraeg, 1991). 

The subject of near wake flows and free 
shear layers in relation to multiple conventional 
anti-armor projectiles in flight is briefly 
discussed by Herrin (1995). In his cited 
example, several of the fired missile projectiles 
are immersed in the near wake of the leading 
warhead (Hohler, 1990). In the context of 
aerodynamics, Herrin emphasizes the 
importance of understanding the wake flow 
physics in order to accurately define the flow 
conditions for the trailing bodies. From these 
few examples, it is evident that wake flows 
found in a variety of practical situations can 
give rise to a diverse range of complex fluid 
phenomenon ranging from flow unsteadiness, 
mixing, separation and flow reversal, to shear 
layer transition and/or reattachment. 

Wake flows at supersonic and hypersonic 
speeds have been studied extensively both 
experimentally and theoretically since the 
1950's. An in depth review of early work and 
discussions of wake characteristics is provided 
by Berger, 1971. In general, hypersonic near 
wake flows are characterized by several 
features which are schematically identified in 
Fig. 2. The inner or near wake flow consists of 
a recirculating fluid that is separated from the 
outer wake flow by a viscous shear layer that 
originates from the forebody boundary layer. A 
weak separation or "lip" shock may form as 
this boundary layer separates from the surface 
(Hama, 1966). As the opposing shear layers 
converge towards the rear stagnation point (in 
the absence of a solid surface), the flow is 
turned through a series of recompression waves 
which eventually coalesce into a recompression 
shock. Flow outside of the shear layer dividing 


2 



streamline continues downstream through a 
neck region; flow inside of the dividing 
streamline is turned back towards the base 
due to recompression. The addition of a solid 
surface will change the flow; but the salient 
features remain similar (Wells, 1990). 

In 1992, a subgroup of an AGAED fluid 
dynamics panel (Working Group 18) was 
formed to address issues concerned with and 
related to blunt body wake flowfields. One of 
the stated experimental objectives of the 
ensuing activity was to examine real gas effects 
in the presence of a strongly expanded wake 
flow behind a 70° blunted cone. Non-intrusive 
flow diagnostics such as Planar Laser Induced 
Fluorescence (PLIF) or Laser Induced 
Fluorescence (LIF) of Nitrous Oxide (NO) were 
to be used to measure rotational temperature 
in the near wake. Surface measurements along 
the model forebody, base, and support sting 
were planned to help quantify blunt-body shear 
layer separation, turning angle, wake 
establishment/ unsteadiness, and wake closure 
in the presence of flows exhibiting real gas 
behavior. The experiments also were to help 
assess and quantify the performance of wind 
tunnels capable of producing real gas flows. 
The tests were to be conducted at a common 
test condition in several of the world's premiere 
hypervelocity shock tunnels: the piston-driven 
shock tunnel at DLR-Gottingen (HEG), the 
Large Energy National Tunnel (LENS) at 
Cubric/Calspan, and the 42-in, shock tunnel at 
NASA Ames (deactivated prior to tests). To 
date, non-intrusive measurements have been 
attempted at HEG but with very limited 
success (Rosenhauer, 1994). Fortunately, 
surface measurements made in these facilities 
have met with greater success and several 
continuum high-enthalpy experimental and 
computational studies have resulted (Kastell 
et. ah, 1994, 1995; Holden, 1994; Holden 

et. al., 1995; Gochberg et ah, 1996). 

It was recognized within the AGARD 18 
activity that along with real gas effects, shear 
layer transition to turbulence would have a first 
order influence on the near wake flow field, and 
in particular, the heating associated with 
reattachment. Furthermore, for a given test 
gas, the separation of these potential viscous 
effects from real gas effects would be difficult, if 
not impossible, to achieve from shock tunnel 
tests alone. A brief review of the relevant 
experimental literature on the subject of shear 


layer reattachment indicated that the peak 
heat transfer associated with an impinging 
transitional shear layer could be much higher 
than when the separated flow is entirely 
laminar or turbulent (Baker, 1966). 
Unfortunately, most blunt-body computational 
studies in the literature do not address the 
issue of shear layer transition but rather focus 
on forebody transition from laminar to 
turbulent flow prior to separation (Mitcheltree, 
1995; Laurien, 1995). In an attempt to 
experimentally identify and separate wake 
viscous effects and wake flow establishment 
factors from potential real gas effects observed 
by Kastell and Horvath (1994), a set of tests 
(on the same model) was conducted at NASA 
Langley at Mach 10 and 6 in two conventional 
low enthalpy blowdown wind tunnels to 
compliment the high enthalpy, hypervelocity 
Mach 10 tests conducted at HEG and LENS. 

The purpose of the conventional 
hypersonic tests at Langley was to obtain both 
flow field and surface measurements for a fully 
established wake flow in a well characterized 
free-stream. The objective of this paper is to 
present results from wake flow field surveys 
conducted in the NASA LaRC 20-Inch Mach 6 
wind tunnel and to determine the state of the 
shear layer and to quantify its position and 
growth. Near wake surface heat transfer 
measurements are presented to evaluate 
viscous shear layer effects on reattachment 
heating. Measured heating distributions are 
compared to laminar Navier Stokes predictions. 

Apparatus and Test 

Model Description 

A photograph and dimensioned sketch of 
the blunt cone model are shown in Figs. 3 and 
4, respectively. The 6-in. diameter 70 deg 
blunted cone forebody was fabricated from 
Chromel rather than stainless steel to avoid 
the potential effects of extraneous EMF's 
resulting from the press-fit installation of 
Chromel-constantan coaxial surface 

thermocouples into a stainless steel model 
(Kidd, 1994; Wieting, 1987). The flat model 
base and 1.5-in. diameter cylindrical support 
sting were machined from conventional steel 
and slotted to accept instrumented ceramic 
(Macor) inserts. Quality assurance 
measurements on the model were made to 
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verify the accuracy of the surface geometry 
(±0.002 -in from nominal) and to precisely 
locate the individual sensors. The forebody 
was instrumented with 40 Chromel-constantan 
coaxial surface thermocouples along a single 
ray, while the base and sting ceramic inserts 
were instrumented with 65 thin film resistance 
gages as described by Miller (1981). Coaxial 
thermocouples were selected for use on the 
forebody because of their small size, fast 
response, and durability from particulate 
damage anticipated during the shock tunnel 
tests at HEG and NASA Ames. 

For some of the tests, aluminum oxide 
grit of approximately 0.381 mm (0.015-in.) 
diameter was dispersed on the forebody 
stagnation region to examine the sensitivity of 
reattachment heating to surface roughness, 
Fig. 5. The ratio of nominal roughness 
diameter to computed boundary layer thickness 
was approximately 3 near the stagnation point 
and decreased to 1.2 at S/R n = 0.9 (end of 
transition grit). 

The uncooled model and support sting 
were attached to the facility support barrel 
approximately 2 model diameters downstream 
in an attempt to maintain setup consistency 
between all tests (HEG, LENS, and LaRC) 
performed on this model. Model angle of 
incidence was fixed at 0 deg. The entire 
assembly were injected to the tunnel centerline 
from a retracted position. 

Facility Description _ 

The model was tested in the 20-Inch 
Mach 6 Air Tunnel at the NASA Langley 
Research Center. A detailed description of this 
facility and related instrumentation is 
presented by Miller (1992). The 20-Inch Mach 
6 is a blowdown facility which uses dried, 
heated, and filtered air as the test gas. Typical 
operating conditions for the tunnel are 
stagnation pressures ranging from 30 to 500 
psia and stagnation temperatures from 76 to 
1000 °R, yielding free stream unit Reynolds 
from 0.5 to 8xl06 6 /ft. It has a closed 20-in by 
20-in. test section with a contoured two- 
dimensional nozzle which provides a nominal 
free stream Mach number ranging from 5.8 to 
6.1. A hydraulically operated model injection 
mechanism can inject the model into the flow in 
0.5 seconds. A photograph of the test section, 
model, and support system arrangement are 


shown, Fig. 6. The maximum run time for this 
facility is approximately 5 minutes; typical run 
times for the separate wake surveys and heat 
transfer measurements for this test series were 
180 and 5 seconds, respectively. 

Mass flow and total temperature 
fluctuations of 2-3 percent and 1-1.5 percent, 
respectively, were measured previously by 
Stainback and Kubendran (1994) in the Mach 
6 test section utilizing a dual wire constant 
temperature anemometer. 

Instrumentation. Data Reduction, and 
Uncertainty 

The reservoir pressure pt ? i was 
measured with two silicon sensors having a full 
scale rating of 500 psia or 150 psia, depending 
on the operating condition of the tunnel. The 
reservoir temperature, Tt,l, was measured 
with two iron-constantan thermocouples 
inserted through the wall of the settling 
chamber. Test-section wall static and pitot 
pressures were monitored and compared to 
tunnel empty conditions to assess model 
blockage effects. Differences in pitot pressure of 
less than 0.5 percent were measured and it 
was concluded that significant blockage did not 
exist. The ratio of projected model frontal area- 
to- tunnel cross sectional area for the present 
test was 0.071. A 16-bit analog-to-digital 
acquisition system acquired the data on all 
channels at a rate of 50 samples per second 
except where noted. 

Near wake pitot pressure was measured 
intrusively using a miniaturized (0.013-in, 
diameter tubing flattened on the probe tip) 
piezoresistive pressure transducer mounted to 
a wedge shaped water-cooled strut (Ashby, 
1988) , as shown in Fig. 7. The close proximity 
of the transducer to the probe tip 
(approximately 3-in.) and a controlled thermal 
environment provided by the water-cooled 
housing jacket minimized pressure settling (lag) 
time and maintained a constant temperature 
at the transducer head during the traverse 
through the wake. Prior to a run, a sealed 
chamber placed over the assembly allowed the 
transducer to be calibrated in-situ. This daily 
in-situ calibration consisted of applying 1 3 
known pressures selected to span the range of 
expected pressures and then applying a second 
order curve fit to the results. The defined 
pressure-voltage relationship was found to be 
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essentially linear. Pitot pressure was 
measured using a transducer rated for a 
maximum pressure of 20 psia to accommodate 
any unanticipated high pressures associated 
with the shear layer recompression process or 
tunnel unstart loads. During the traverse, the 
probe was stepped in increments of 0.003-in. 
with a dwell point at each position selected to 
allow the local pressure to settle. The probe 
position was determined from a known 
reference point on the model and the traverse 
stations (for both hot-wire and pitot surveys) 
preprogrammed prior to a run. While no 
attempt was made to correct for misalignments 
due to wind on loading during a run (model and 
survey support systems are independent), 
schlieren observations suggested little relative 
movement. 

Schlieren videos (30 frames/sec) and still 
photographs (spark duration of 5nsec) were 
made of each run for visualization of the flow 
field. When necessary, digital enhancements of 
the schlieren images were performed using 
commercial software to bring out details of the 
near wake flow field. 

A constant temperature hot wire 
anemometer was utilized to provide fluctuating 
voltage measurements and spectra to assess 
the state of the shear layer. The locations of 
the hot wire surveys were selected to 
correspond to the pitot measurement stations. 
Time histories of the instantaneous voltage 
fluctuations provided qualitative information 
regarding the distribution of r.m.s. energy. The 
constant temperature anemometer consisted of 
a 5 pm diameter Platinum-Rhodium wire 
operated at an overheat ratio of approximately 
1 and sampled at 500 kHz. A more complete 
description of the acquisition system and 
analysis techniques can be found in Spina 
(1994). 

Forebody surface temperature 

measurements were obtained using 
commercially available 0.031-in. diameter 
Chromel-constantan coaxial surface 

thermocouples that were press fit into the 
model wall. The model contoured junction was 
formed at the sensing surface by blending 
(lightly sanding) the two materials together. 
Several thermocouples had a second junction 
that was used to monitor the backside wall 
temperature in order to assess conduction 
effects. A general description of the coaxial 


thermocouple and its use is outlined by Kidd 
et. al. (1994). 

The more sensitive thin film resistance 
gage was used to measure surface temperature 
in the base and wake regions. Standard 
mechanical deposition techniques developed at 
LaRC (Miller, 1981) were used to fabricate the 
0.030-in. by 0.040-in. platinum sensing 
element. Surface temperatures were integrated 
over time to determine the local heat transfer 
rate using the computer code developed by 
Hollis (1995). Both analytical (Cook, 1970; 
Kendall-Dixon, 1967) and numerical finite- 
volume heat transfer models are incorporated 
into this code. The analytical solutions are 
derived from one-dimensional, semi-infinite 
solid heat conduction theory with the 
assumption of constant substrate (model) 
thermal properties. When using the analytical 
option the inferred heating rates are empirically 
corrected for the effects of variable model 
thermal properties. The finite-volume 
technique, which was used for the present 
tests, directly accounts for the variable model 
thermal properties (results from recent 
laboratory tests to determine substrate 
thermal properties have been incorporated into 
the code; Hollis, 1995). This option also 
removes the restriction of a semi-infinite wall 
boundary condition. This was an important 
consideration for this study as the model was 
designed and instrumented primarily for 
impulse facility testing where test times are on 
the order of milliseconds and the semi-infinite 
wall assumption is valid. In contrast, test 
times in the conventional hypersonic wind 
tunnels at LaRC were several orders of 
magnitude longer and led to invalidation of this 
assumption. For the present study, the 
uncertainty associated with variable wall 
thermal properties is believed to be minimal, 
particularly in the wake region where surface 
temperature increases of 30 °F or less were 
measured. 

A typical temperature and Ch time 
history for a thin-film sensor located along the 
support sting in the wake region is shown in 
Fig. 8. Unlike shock tunnel data where wake 
flow establishment times must be determined, 
the local Stanton number remains essentially 
constant upon reaching the tunnel centerline. 

Measured values of Pt 5 i and Tt,l are 
believed to be accurate to within ±2 percent. 
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Based on manufacturer specifications, values of 
the near wake pitot pressure are estimated to 
be accurate to within ±5 percent. Based on 
previous experience, the accuracy of the heat 
transfer measurements is believed to be better 
than ±10 percent. Repeatability for the heat 
transfer measurements was found to be 
generally better than ±1 percent. 

Computational method 

The NSHYP code, used in the present 
study, was developed at DLR for the 
computation of supersonic and hypersonic flows 
of a perfect gas or a chemically reacting mixture 
of perfect gases (see Brenner et al., 1993, 
Brenner and Prinz, 1992, Riedelbauch and 
Brenner, 1990). Many hypersonic flow fields of 
interest permit the use of the thin-layer Navier- 
Stokes equations. However, due to the large 
recirculation region which develops in the wake 
of the blunted cone, the basic equations under 
consideration here are the full unsteady Navier- 
Stokes equations for a body-fitted coordinate 
system (x: streamwise, h: circumferential, z: 

wall normal coordinate), 

dJJ dE dF dG 1 

dt d% dt] d£ Re 

» ,r 

where 

U - J-'(p,pu,pv,pw,ej 

The inviscid and viscous fluxes are denoted by 
E y F>G and E ,F ,G respectively. The 
transformed vector of the conservative variables 

is given by U ; here J 1 represents the jacobian 
of the transformation matrix. The equations 
have been nondimensionalized using the free 
stream density, velocity and viscosity and the 
radius of the blunted cone. For perfect gas 
computations the viscosity coefficient m and the 
coefficient of thermal conductivity k are 
calculated from the Sutherland law. The 

Prandtl number is assumed to be constant. 
The thermodynamic properties pressure and 
temperature are calculated from the 

conservative variables using the thermal and 
caloric state equation of perfect gas with the 


ratio of specific heats for diatomic gases (g 

= 1.4). 

The time discretization is fully implicit in 
order to avoid stiffness problems introduced by 
the small mesh increments in the boundary 
layer. NSHYP was primarily designed for the 
computation of steady flow fields. Therefore, 
the time integration was chosen to be of first 
order accuracy as only the converged steady- 
state solitions are of interest. The spatial 
derivatives of the inviscid fluxes are 
approximated by a second order upwind Total 
Variation Diminishing (TVD) formulation 
according to Yee and Harten (1987), and the 
spatial derivatives of the viscous fluxes are 
discretized by second order central differences. 
The TVD flux-difference splitting algorithm 
involves the solution of locally one-dimensional 
Riemann problems at the cell interfaces. Here 
the approximate Riemann solver of Roe (1981) 
is applied. Introducing a line Gauss-Seidel 
relaxation results in a linear system of block- 
tridiagonal matrices which is solved with the 
Richtmyer algorithm. 

The physical domain considered here is 
bounded by the body, the inflow and outflow 
boundary and the line of symmetry. On the 
body no-slip conditions apply. The temperature 
of the isothermal wall is set to T wa n = 300K. 
The momentum equation in wall-normal 
direction is approximated by dr/dn = 0 at the 
wall. At the inflow boundary, a homogeneous 
supersonic flow is assumed. The free-stream 
conditions are given in Table 1 . At the outflow 
boundary the conservative variables are 
extrapolated from the integration domain by 
assuming that their slope in the x-direction is 
constant, i.e. d 2 U/dx 2 = 0. In order to compute 
axisymmetric flows with the present 3D code, 
appropriate symmetry conditions were 
employed to evaluate the fluxes in 
circumferential direction. At the line of 
symmetry, which for axisymmetric flows 
represents a singularity, the variables were 
determined by means of a cubic extrapolation. 

Grid description and sensitivity studies 
will be deferred to the discussion of the 
comparisons of measurement with prediction. 

Test conditions 

Tests were performed at Mach 6 in air for 
nominal reservoir pressures and temperatures 
shown in Table 1. The nominal free stream 
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Reynolds numbers based on model diameter for 
these conditions were 0.5xl0 6 , 2xl0 6 , and 
4xl0 6 . A majority of the tests were conducted 
at Reoo ? d = 2xl0 6 (lower dynamic pressure) to 
minimize the chances of breaking the hot wire 
while maintaining high enough wake density 
gradients to visualize the separated shear layer 
via schlieren. 

Results and Discussion 

Approach 

Initial tests were conducted on the 
blunted cone in the NASA LaRC 31 -Inch Mach 
10 Tunnel prior to the planned high enthalpy 
tests at HEG or LENS. These tests revealed 
the sensitivity of wake shear layer 
reattachment heating to Reynolds number, 
Fig. 9. The magnitude of the heating peak 
associated with reattachment is in the range of 
16-18 percent of the measured forebody 
stagnation point heating These wake peak 
heating levels were also measured by Dye 
(1993) on a similar 70 deg sphere cone tested 
in the same facility. However, subsequent 
shock tunnel data (Holden, 1995) obtained on 
the present blunted cone configuration and 
corresponding laminar predictions from Moss 
(1995) imply that peak heating associated with 
laminar shear reattachment is typically 5-6 
percent of the forebody stagnation point 
heating. One explanation for the differences 
observed between tests in the 31-Inch Mach 10 
tunnel and LENS is that the relatively higher 
LaRC wake heating was produced by a non- 
laminar reattachment process. 

The tests in the NASA LaRC 20-Inch 
Mach 6 Tunnel were initiated to provide both 
flow field and surface measurements to 
complement the Mach 10 results. The original 
objective of the experiment was to spatially 
locate the shear layer and define its growth via 
optical and intrusive measurements. The state 
(laminar/transitionahturbulent) of the shear 
layer and its influence on reattachment heating 
was to be determined using hot wire 
anemometry and surface heat transfer 
measurements, respectively. A condition was 
sought corresponding to a laminar shear layer. 
If a laminar condition was identified, shear 
layer transition would then be promoted using 
localized roughness placed on the forebody. 


Based on similar magnitudes in the 
measured ratio of sting peak heating to 
forebody stagnation point heating from the 
present Mach 6 tests, a laminar condition was 
not identified. Because of this situation, CFD 
predictions were used to characterize laminar 
reattachment heating levels. The experimental 
measurements were used to establish the 
existence of a transitional/turbulent shear layer 
and to provide heating measurements for this 
type of reattachment process. 

The data presented in the subsequent 
sections are organized and presented in the 
following manner: (1) Shear layer thickness vs. 
axial position are presented over the range of 
Reynolds numbers. (2) Surface heating 
distributions along the model forebody, base, 
and support sting are presented in the form of 
a normalized Stanton number, Ch/Ch,stag> 
where Ch,stag corresponds to the measured 
stagnation point heating on the model forebody 
unless otherwise noted. (3) Hot wire 

fluctuating voltage time histories and spectra 
are highlighted to support the hypothesis that 
a transitional/turbulent shear layer exists prior 
to reattachment, and (4) Comparisons of 
measured heating with prediction along the 
forebody, base, and sting are made in the form 
of an unnormalized Stanton number, Ch- 

Pitot pressure surveys 

Near wake pitot pressure profiles were 
used to locate the mean shear layer 
boundaries. Based on schlieren photographs, 
several wake axial survey stations of interest 
were identified, Fig. 10. An example of a 
typical pitot pressure profile along a fixed axial 
station is shown, Fig. 11a. A schlieren 
photograph, Fig. lib, illustrates the 
corresponding general flow features. As 
expected, the pressure decreases as the probe 
descends through the expansion fan. The small 
jump in the pressure near the trailing edge of 
the expansion fan corresponds to the lip shock 
discussed earlier (Hama, 1966). In this study, 
the boundary of the free shear layer is 
determined by the local slope change in the 
pitot pressure profile; visually, the shear layer 
can be seen in the schlieren image. 

Shear layer boundaries inferred from a 
series of pressure profiles obtained along the 
wake axis are superimposed on a schlieren 
image obtained from a separate run, Fig. 12. 
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The close proximity of the mean shear layer 
boundaries determined via intrusive methods 
to that obtained optically, suggests that 
interference effects which could be expected 
from the survey strut and probe were not large. 
At low pressure conditions, where the density 
gradients were not sufficiently large enough to 
visualize the shear layer, these intrusive 
measurements were essential to properly 
locating the hot wire. 

Shear layer thickness as a function of 
axial distance is generally linear as shown in 
Fig. 13. The measured shear layer thickness 
increases with decreasing Reynolds number. 

Surfa ce heating 

Surface heating distributions obtained 
with and without forebody transition grit are 
shown over a range of Reynolds numbers, 
Fig. 14a-c. As expected for the smooth 
forebody, the heating in the stagnation region 
rapidly diminishes as the flow accelerates 
around the spherical nose. As the boundary 
layer expands and thins near the forebody 
shoulder the heating increases and then quickly 
falls off as the boundary layer separates from 
the base. On the model base low levels of 
heating are measured and further downstream 
a local peak occurs that is associated with the 
impingement of the shear layer on the sting. 

The influence of roughness (via the 
transition grit) on the forebody and near wake 
heating is particularly evident at the two 
highest Reynolds numbers, Fig. 14b-c. From 
these data it is concluded that at the two 
highest Reynolds numbers (Reoc,d = 2xl0 6 and 
Reoc^d = 4xl0 6 , the forebody boundary layer has 
transitioned to turbulent flow. The associated 
forebody heating has been augmented by up to 
a factor of two. This observation is consistent 
with flight calculations reported by Mitcheltree 
(1995) for laminar/turbulent forebody convective 
heating to a 70 deg blunted cone similar to the 
model tested. At Reoc,d = 0.5xl0 6 , Fig. 14a, 
the roughness was insufficient to promote a 
turbulent boundary layer and the forebody 
heating returned to nearly laminar levels as the 
flow expanded around the model shoulder. 

The effect of Reynolds number on heating 
is shown, Fig. 15 a-b, for the smooth and rough 
models, respectively. In both situations, the 
heating peak associated with the near wake 
reattachment process increases in magnitude 


and moves towards the base with increasing 
Reynolds number. This behavior is indicative 
of a transitional/turbulent shear layer (Holden, 
1995) (Berger, 1971). 

Fig. 16 a-b presents the forebody heating 
distributions as measured on the rough model 
up to the point of separation on the model 
shoulder. The first data set, Fig. 16a, has 
been normalized by the measured stagnation 
point heating on the rough model. The second 
heating distribution, Fig. 16b, is the identical 
data normalized to the stagnation point 
heating as measured on the smooth model. 
The first format is selected to more clearly 
discuss trends, while the latter is selected to 
illustrate the apparent stagnation point 
heating augmentation due to the roughness 
located in the stagnation region. Fig. 16a 
indicates that the forebody boundary layer has 
transitioned to turbulence at S/R n =0.4 with an 
associated heating maximum of approximately 
35 percent above that measured at the 
stagnation point on the roughened model. The 
tripped heating distributions suggest a 
relaminarization of the boundary layer along 
the cone as it approaches the model shoulder at 
S/R n = 2.1. The tripped flow found at the 
higher Reynolds numbers appears to have 
diminished the small local heating peak 
observed at the shoulder (S/R n = 2.1) for 
laminar conditions, Fig. 15a. Noting this 
same trend in his turbulent flight calculations, 
Mitcheltree (1995) questioned the creditability 
of his turbulence model in this region of strong 
expansion. The present perfect gas data 
implies that the calculation of Mitcheltree in 
the shoulder region might be valid. 

Fig. 16b normalizes the rough wall 
heating measurements to the smooth wall 
stagnation point heating. At the stagnation 
point, a heating augmentation of 17 percent 
above laminar levels at Reoo ? d = 4xl0 6 is noted. 
Holden (1983), reports several instances from 
wind tunnel and flight tests where boundary 
layer transition occurred in or just downstream 
of the stagnation region. In these examples 
cited by Holden, the stagnation point heating 
was found to be significantly larger than those 
predicted on the basis of laminar theory alone. 
He performed several studies to examine 
potential mechanisms which may have resulted 
in the enhanced stagnation heating. One 
study, involved trip rings placed 10 (or more) 
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boundary layer thicknesses downstream of the 
stagnation point. No significant heating 
increase at the stagnation point was measured. 
The roughness elements in the present study 
were placed immediately adjacent to the 
stagnation point (k/d = 3 at S/R n = 0.05) and it 
is apparent that some form of influence is 
present for this situation. 

Measured near wake heating 
distributions along the model base and sting 
over the range of Reynolds numbers with and 
without forebody transition grit are presented, 
Fig. 17a-b. The heating peak associated with 
the reattachment process is evident and occurs 
1 to 1.3 model diameters (S/R n = 7.8 to 8.8) 
downstream of the model base. Reattachment 
heating associated with the smooth model, 
Fig. 17b, indicate that the magnitude of the 
peak was approximately 21-29 percent of the 
measured forebody stagnation point suggesting 
a non laminar reattachment process. Tripping 
the forebody boundary layer tended to reduce 
the peak heating associated with 
reattachment, Fig. 17a. Similar to LaRC Mach 
10 data on the same model, Fig. 9 and results 
presented by Hollis, (1995), an increase in 
Reynolds number produced an increase in 
magnitude and forward movement (towards the 
base) of the wake heating maximum. 

Holden (1995) speculated that differences 
in measured wake heating levels from a broad 
range of AGARD 18 blunt-body tests could 
result from differences in the turbulent 
structure of the shear layer. As discussed 
earlier, Berger (1971) discussed several wake 
flow models that suggest the size of a laminar 
wake recirculation region should increase with 
increasing Reynolds number. The present 
Mach 6 wake heating trends imply that shear 
layer transition has occurred prior to 
reattachment/recompression. The heating on 
the model base remains below 6 percent of the 
stagnation level but is quite sensitive to 
Reynolds number. Complex recirculation 
regions have been noted in the present 
computations which may result in the local 
heating peaks that are measured on the model 
base. These results will be discussed in the 
Comparisons to Measurement section. 

The influence of forebody transition grit 
on wake heating is summarized, Fig. 18. At 
Reoo ? d = 0.5xl0 6 there is no influence of the 
forebody roughness on the reattachment 


heating. The forebody heating distribution at 
this condition indicates boundary layer 
relaminarization prior to separation at the 
blunted cone shoulder. Based on the level of 
heating at reattachment, it is believed that the 
shear layer is transitional at the point of 
impingement on the sting. At Reoc,d = 2xl0 6 
and 4xl0 6 the peak heating measured on the 
sting was reduced by approximately 20 percent 
and 24 percent, respectively, when the forebody 
transition grit was applied. The spatial 
location of the heating maximum remained 
unchanged. Several repeat wind tunnel runs 
were to verify this trend: the experimental 
heating distributions were repeatable to within 
1.1 percent. The forebody heating distribution 
at this condition indicates boundary layer 
transition has occurred prior to separation at 
the blunted cone shoulder. It is believed that 
at reattachment the free shear is turbulent. 

These trends can be compared with 
supersonic heat transfer measurements made 
by Naysmith (1961) and later by Baker (1966) 
for reattaching flow over rearward facing steps. 
In these early studies, it was concluded that 
when transition from laminar to turbulent flow 
occurred along the separated shear layer, the 
peak heat transfer in the reattachment zone 
was much higher than when the shear layer 
was completely laminar or turbulent. The 
implication to the present results suggest that 
at a high enough Reynolds number, without 
the forebody transition grit, the instability 
mechanisms present (ie. tunnel noise, surface 
roughness etc.) were insufficient to produce a 
fully developed turbulent shear layer. The 
addition of surface roughness on the forebody 
may have been sufficient to promote a more 
developed turbulent shear layer which resulted 
in a heating reduction. 

Hot wire surveys 

To further clarify the state of the shear 
layer, an uncalibrated hot wire was traversed 
at stations corresponding to the pitot surveys. 
The intention was to try to detect laminar or 
turbulent signals in the shear layer and the 
data is strictly qualitative. Traverses were 
obtained at Reoo t d = 0.5xl0 6 and Reoo ? d 
= 2xl0 6 . A sequence of fluctuating voltage time 
histories obtained at Reoc,d = 2xl0 6 and x/d 
= 0.333 for specific points within the shear 
layer is shown in Fig. 19a-d. Outside of the 
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shear layer Fig. 19a, the fluctuations are small. 
As the wire enters the shear layer, Fig. 19b, 
the signal becomes intermittent with occasional 
voltage drops which indicate times when the 
flow is decelerated. This behavior is consistent 
with turbulent signals when entering a shear 
layer. Toward the middle of the shear layer, 
Fig. 19c, the flow fluctuations increase in 
amplitude. As the wire departs the shear 
layer, Fig. 19d, headed toward the recirculating 
region, a pronounced intermittancy returns 
except now the voltage spikes are positive. 
Figure 20 shows the power spectra associated 
with trace (a) and (c) from the previous figure. 
The spikes in the data above 30kHz are from 
strain gauging and electronic noise and do not 
contribute substantially to the r.m.s levels. The 
spectra indicates that there is a spread of 
energy over a large frequency range. 
Furthermore, the energy in the middle of the 
shear layer is considerably larger than that 
outside the shear layer. These behaviors are 
characteristic of turbulent signals. Similar 
surveys were taken at x/d = 0.042, and 0.500. 
The data taken at x/d = 0.042 suggest that the 
flow is at least transitional at this station. All 
surveys indicate that the r.m.s voltage 
normalized by the mean voltage increases as 
the center of the shear layer is approached. 
Traveling downstream along the shear layer, 
the peak r.m.s levels increase quickly between 
x/d = 0.042 and x/d = 0.333 and then begin to 
level off between x/d = 0.333 and x/d = 0.500. 
This may be indicative of a shear layer 
completing the transition process. Hot wire 
data obtained at Reoo ? d = 0.5xl0 6 display 
similar trends but with much lower r.m.s 
levels. 

Computational Predictions 

Grid sensitivity 

The computational results from the 
present study indicate the necessity of grid 
refinement studies to examine the sensitivity of 
the wake flow field and associated surface 
heating to grid spacing; particularly regarding 
blunt body base flow calculations. The topology 
of the nominal single block grid used in the 
present computations is shown, Fig. 21. In 
order to investigate the dependence of the 
solution on grid resolution, three grids 
containing 386x151, 193x76, and 97x38 points 


in wall tangential and wall normal direction, 
respectively, were used. The finest grid was 
first generated using the elliptic grid generator 
MEGACADS (Ronzheimer et. al. 1994) For this 
grid, the distance of the first grid line to the 
wall normalized to the radius of the forebody 

was 6xl0" 6 . The remaining two grids were 
obtained by removing every second grid point in 
each direction. 

The effect of the grid refinement on 
surface heating for Reoo,d = 2xl0 6 are 

summarized in Fig. 22 a-b. All solutions 
obtained with the three grids were steady. As 
the comparisons of experimental data with 
prediction in this paper are limited to surface 
heat transfer, the discussion of grid convergence 
is restricted to the effect on heating 
distributions. The computed surface heat 
transfer along the forebody, Fig. 22a., showed 
little influence to grid refinement except at the 
stagnation point and on the model shoulder. 
The predicted heating at these two locations, 
however, was reduced with each successive 
refinement. For example, solutions obtained at 
these locations with the medium and fine grid 
agree to within 2 percent at the stagnation 
point and 8 percent at the model shoulder. 

The predicted surface heating along the 
sting obtained with the three grid sizes, 
Fig. 22b, indicated that the heating peak 
associated with shear layer reattachment 
increased and moved slightly foreword towards 
the model base. Solutions obtained with the 
medium and fine grid show the magnitude of 
this heating peak remained essentially 
constant while the location was displaced 3 
percent. Particle traces associated with the 
three grids, Fig. 23a-c, reveal the complexity of 
the predicted flow field found at the model base 
sting juncture. The particle traces shown were 
obtained by integration of the steady velocity 
field. Refinement of the grid appeared to 
resolve an additional fourth vortex within the 
base recirculation region (Fig. 23b-c). The 
reattachment point associated with the small 
vortex found just behind the model shoulder 
remained spatially fixed for all three grids. 

In addition to the effects of grid 
refinement the sensitivity of surface heating to 
numerical damping was also investigated. The 
upwind TVD scheme which was used to 
perform the present computations requires an 
entropy correction to produce a numerically 


10 



stable solution. In this study, this form of 
artificial viscosity was varied from the 
maximum value permitted (see Yee, 1987) to 
the minimum value which still allowed a stable 
solution to be computed across the bow shock 
wave. The results from this investigation (not 
shown) indicated that the stagnation point 
heating obtained with the fine grid increased by 
6 percent as the entropy correction parameter 
was reduced to its minimum value. At the 
same time, the localized heating peak found at 
the model shoulder decreased by approximately 
6 percent. These percentages do not truly 
represent an expected uncertainty as the large 
numerical damping utilized in this parametric 
study are not routinely employed. The 
computed wake heating and base flow features 
were not affected as the artificial viscosity 
levels were varied. 

Based on the present grid sensitivity 
study it has been concluded that accurate 
solutions have been achieved. The small 
differences in heating predicted between the 
medium and fine grid at the stagnation point 
and model shoulder are comparable to 
conclusions from a blunt body grid sensitivity 
study performed by Moss et. al. (1993). From 
an engineering perspective, the computed 
heating distribution obtained with the fine grid 
represents a conservative estimation. 
Extrapolation to finer meshes would yield lower 
predicted heating at the forebody stagnation 
point and shoulder. The magnitude of the 
shear layer reattachment heating peak on the 
sting did not change between solutions 
obtained with the medium and fine grids. 

Comparisons to Measurement 

Comparisons of predictions obtained with 
the fine mesh (386 x 151 grid) to 
measurements along the model forebody and 
base-support sting are shown in Figs. 24a-b, 
25a-b, 26a-b, for Reoc,d = 0.5xl0 6 , 2xl0 6 , and 
4xl0 6 respectively. An unnormalized Stanton 
number was selected for comparative purposes 
as normalization can often mask differences 
between the distributions. On the forebody, 
Figs 24a, 25a, and 26a, the best agreement 
(less than 1.5 percent difference) was found 
along the cone section for Reoc,d = 0.5xl0 € and 
Reoo,d = 2xl0 6 . The largest disparity between 
measurement and laminar prediction on the 
cone section was 9 percent at Recced = 4xl0 6 . 


At this higher Reynolds number it is possible 
the forebody boundary layer was approaching a 
transitional state. Although the measurements 
shown in Figs. 24a, 25a, and 26a were 

obtained without the use of transition grit, the 
particulate induced roughness on the model 
from the prior HEG shock tunnel tests may 
have been sufficient to promote a forebody 
boundary layer instability at Reoo ? d = 4xl0 6 . 

Spatial resolution of the heat transfer 
gages in the vicinity of the model shoulder was 
inadequate to resolve the local heat peak 
observed here computationally. The 

uncertainty of the measured heat transfer is 
greatest at the shoulder due to thermal 
conduction effects not accounted for in the data 
reduction process. The asymptotic behavior of 
the predicted heating as it nears the stagnation 
point of the spherical nose suggests that the 
finer, 386 x 151 grid was required to resolve 
the boundary layer and associated gradients. 
The use of this grid reduced the difference 
between stagnation point measurement and 
prediction to 2.3 percent at Recc,d = 4xl0 6 . 

The laminar calculations in the near 
wake indicated that the heating peak 
associated with reattachment was spatially 
located downstream of that inferred from 
measurement and was approximately 50 
percent below the measured value in 
magnitude, Figs 24b, 25b, and 26b. This is 
consistent with the conclusions drawn earlier in 
this study. That is, for a transitional/turbulent 
reattaching shear layer, the size of the 
recirculation region is expected to diminish and 
the magnitude of the heating peak to increase 
relative to the laminar counterpart. The 
computed heating peak found downstream on 
the sting varied with Reynolds number and 
was 10-14 percent of the predicted forebody 
stagnation point heating. The corresponding 
experimental results varied from 21-29 percent. 

Generally, the predicted base heating 
was twice that measured. It is assumed the 
local peak found on the base (S/R n = 2.3) is the 
result of flow stagnation associated with the 
small vortex near the shoulder observed 
computationally, Fig. 23a-c. The predicted 
base heating associated with this flow 
structure was sensitive to Reynolds number 
and varied from 4-12 percent of the predicted 
forebody stagnation point heating. Validation 
of the complex recirculation patterns predicted 
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near the base-sting juncture were not possible 
as surface instrumentation was not placed in 
this location. 


Concluding Remarks 

The reattachment process of an 
axisymmetric free shear layer associated with 
the near wake of a blunt-body was 
experimentally and computationally 

investigated at Mach 6 over a free-stream 
Reynolds number range based on body 
diameter of 0.5xl0 6 to 4xl0 6 . The configuration 
consisted of a 70 deg half angle, spherically 
blunted cone, with a cylindrical after body. 
Surface heating rates were computed from 
temperature-time histories from coaxial surface 
thermocouples on the model forebody and thin 
film resistance gages along the model base and 
cylindrical after body. The general flow 
features (bow shock, wake shear layer, and 
recompression shock) were visually identified by 
schlieren photography. Mean shear layer 
position and growth were determined from 
intrusive pitot pressure surveys. In addition, 
wake surveys with a constant temperature hot- 
wire anemometer were utilized to qualitatively 
characterize the state of the shear layer prior to 
reattachment. Perfect gas numerical 
simulations were provided by a 3-D Navier 
Stokes code (NSHYP). 

Based on the present experimental 
results, the following conclusions can be made. 
Shear layer impingement on the instrumented 
cylindrical after body resulted in a localized 
heating maximum that was 21 to 29 percent of 
the forebody stagnation point heating. Peak 
heating resulting from the reattaching shear 
layer was found to be a factor of 2 higher than 
laminar predictions, suggesting a transitional 
or turbulent shear layer. An increase in 
Reynolds number produced an increase in 
magnitude and forward movement (towards the 
base) of the wake heating maximum on the 
sting; behavior indicative of a turbulent wake 
shear flow. Fluctuating voltage time histories 
and spectra from the hot wire surveys across 
the shear layer support this observation. 

The sensitivity of surface heating to 
forebody roughness was characterized for a 
reattaching transitional shear layer. For 


example, at Reoo,d = 4xl0 6 , when the shear 
layer was transitional, the magnitude of peak 
heating from shear layer impingement was 
reduced by approximately 24 percent when 
transition grit was applied to the forebody. 
The spatial location of the local peak, however, 
remained unchanged. 
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Table 1. Nominal Flow Conditions 


Reoo,d 

0.5 xlO 6 

1 xlO 6 

2 xlO 6 

3.5 xlO 6 

P t ,i (N/ m 2) 

4.05 xlO 5 

8.62 xlO 5 

17.2 xlO 5 

32.7 xlO 5 

T t ,l (°K) 

493.7 

508.8 

505.0 

523.9 

roo (kg/m^) 

1.68 xlO' 2 

3.22 xlO' 2 

6.45 xlO" 2 

11.54 xlO' 2 

Too (°K) 

63.3 

63.4 

62.7 

64.4 

uoo (m/sec) 

931 

947.6 

943.8 

962.2 

Moo 

5.84 

5.93 

5.94 

5.98 
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Fig. 1. Conceptual drawing of an aerobraking vehicle . 



Recirculating ■ Rear 

flow stagnation point 


Fig. 2. Schematic representation of the base flow 
region behind a blunt body in hypersonic flow. 



Fig. 3. Photograph of instrumented 70° blunted 
cone heat transfer model. 


I 8 = 18-ln 



Fig. 4. Dimensional sketch of 70° blunted 
cone model. 



Fig. 5. Transition grit dispersed on blunted 
cone (front view). 



Fig. 6. Blunted cone and pitot sur\>ey system installed 
in LaRC 20-in. Mach 6 test section. 


17 










c h /c h, stag >. Temp(°F) 


t ~ ■ 

\ 


Heat-transfer 

gages 


0.013-in 

O.D. 


Warier cooled 

composite 

jacket houses 

miniature 

piezo-resistive 

pressure 

transducer 


Fig. 7. Pitot survey assembly . 



5 6 7 B 9 

Timefsec) 

ig . 5. Temperature and Stanton number time 
history for a wake thin film gage. 




Fig. 10. Schlieren photograph indicating pitot pressure 
and hot wire axial survey stations. 

A/oo = Re oo d ~ 




(b) Schlieren photograph 


Fig. 9. Effect of Reynolds number on normalized 
sting stanton number heating distribution 
A/oo ~ 10, smooth forebody. 


Fig. II. Typical pitot pressure distribution and general 
flow features observed from blunt body near wake 
flowfield. = 6, Re x j = 2x1 (ft Sta. B 
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Shear layer thickness, in. 



Fig. 12. Comparison of shear layer thickness inferred 
from pitot pressure surveys with schlieren image. 
Woo = 6, R eoo d = 4xl(P 



x/d 


Fig. 13. Effect of Reynolds number on shear layer 
thickness as measured from wake pitot surveys. 
M 00 = 6 



(a) A/oo = 6, Re^ j = 0.5xl(fi (k ~ 0.381 mm) 




s/R n 

(c) = 6, Re^ d = 4x1 (P (k ~ 0.381 mm) 


Fig . 14. Effect of forebody roughness on normalized 
Stanton number heating distribution. 
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(a) Af^ = 6, smooth fore body 



s/R n 

( b ) A/oq = 6, rough forebody (k ~ 0.38 J mm) 


Fig. 15. Effect of Reynolds number on normalized 
Stanton number heating distribution. 




(b) = 6, (Cfr sta g smooth forebody) k ~ 0.38 J mm 

Fig. 16. Effect of Reynolds number on normalized 
forebody Stanton number heating distribution. 
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c h /c h, stag c h /c h, stag 




3.0 4.0 5.0 6.0 7.0 8.0 9.0 10 


S/R n 

(b) M oo = 6, smooth forebody 

Fig. 17. Effect of Reynolds number on normalized 
base/sting Stanton number heating distribution. 



0.20 0.40 0.60 0.80 1.0 1.2 1.4 1.6 


x/d 

Fig. 18. Effect of forebody roughness on normalized 
sting Stanton number heating distribution . 

Afoo = 6 y (k ~ 0.381 mm) 
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Power, 




Fig. 19. Fluctuating voltage time histories during a typical hot wire traverse thru shear layer. 


f L 10- 



10 10 10 10 10 
Frequency, Hz 

Fig 20. Power spectra outside and within the shear layer. 
M oo = 6, Re oo j = 2x1 Cp , x/d = 0.333, 
rough forebody 
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(a) Forebody 
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193x76 
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^ n 


Fig 21. Computational grid for NS code simulation 
of blunt body flow. Mesh size 193 x 76. 


[Sting 


3.0 4.0 5.0 6.0 7.0 8.0 9.0 10 

S/Rn 

(b) Base- sting 

Fig 22. Effect of grid refinement on predicted Stanton 
number heating distribution. 

A/oo = 6, Re^ = 2x1 (P, laminar 
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(a) Fore body 



Fig 24. Comparison of measured Stanton number 
heating distribution with prediction. 

= 6, Re x d = O.SxJO 6 


Fig 23. Particle paths for computational wake flow fie Id. 
= 6, Re^ a = 2xl(P y laminar 
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Fig 25. Comparison of measured Stanton number 
heating distribution with prediction . 

M 0 o = 6 f Re^ - 2xl(P, laminar 


Fig 26. Comparison of measured Stanton number 
heating distribution with prediction. 

A/oo = 6, R eoo d = 4x1 (P 
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